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Abstract 

Mercury's capture into the 3:2 spin-orbit resonance can be explained as a result of its 
chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo 
spatial missions is to accurately measure Mercury's rotation and its obliquity in 
order to obtain constraints on internal structure of the planet. Analytical approaches 
at the first order level using the Cassini state assumptions give the obliquity constant 
or quasi-constant. Which is the obliquity's dynamical behavior deriving from a 
complete spin-orbit motion of Mercury simultaneously integrated with planetary 
interactions? We have used our SONYR model (acronym of Spin-Orbit iV-bodY 
Relativistic model) integrating the spin-orbit iV-body problem applied to the solar 
System (Sun and planets). For lack of current accurate observations or ephemerides 
of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical 
integration, we have built an original method for finding the libration center of 
the spin-orbit system and, as a consequence, for avoiding arbitrary amplitudes in 
librations of the spin-orbit motion as well as in Mercury's obliquity. The method has 
been carried out in two cases: (1) the spin-orbit motion of Mercury in the 2-body 
problem case (Sun-Mercury) where an uniform precession of the Keplerian orbital 
plane is kinematically added at a fixed inclination (S2K case), (2) the spin-orbit 
motion of Mercury in the A-body problem case (Sun and planets) (S n case). 

We find that the remaining amplitude of the oscillations in the S n case is one order 
of magnitude larger than in the S2K case, namely 4 versus 0.4 arcseconds (peak-to- 
peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for 
a difference of 10.8 arcseconds. These theoretical results are in a good agreement 
with recent radar observations but it is not excluded that it should be possible to 
push farther the convergence process by drawing nearer still more precisely to the 
libration center. We note that the dynamically driven spin precession, which occurs 
when the planetary interactions are included, is more complex than the purely 
kinematic case. Nevertheless, in such a iV-body problem, we find that the 3:2 spin- 
orbit resonance is really combined to a synchronism where the spin and orbit poles 
on average precess at the same rate while the orbit inclination and the spin axis 
orientation on average decrease at the same rate. As a consequence and whether it 
would turn out that there exists an irreducible minimum of the oscillation amplitude, 
quasi-periodic oscillations found in Mercury's obliquity should be to geometrically 
understood as librations related to these synchronisms that both follow a Cassini 
state. Whatever the open question on the minimal amplitude in the obliquity's 
oscillations and in spite of the planetary interactions indirectly acting by the solar 
torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium 
state that proceeds from a 2-body Cassini state. 

Key words: Mercury, rotational dynamics, spin-orbit resonance, celestial 
mechanics 
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1 Introduction 



From radar measurements of the rotation of Mercury obtained by Pettengill 
& Dyce (1965), Colombo (1965) showed the non-synchronous rotation of Mer- 
cury and its 3:2 spin-orbit state (the rotational and orbital periods are 56.646 
and 87.969 days respectively). The stability of the 3:2 spin-orbit resonance de- 
pends on a significant orbital eccentricity and the fact that the averaged tidal 
torque trying to slow the planet further is considerably less than the maxi- 
mum averaged torque acting on the permanent axial asymmetry that causes 
the planet spin to librate about the resonant 3:2 value. Moreover, Correia & 
Laskar (2004) have shown more recently that Mercury's capture into the 3:2 
resonance can be explained as a result of its chaotic orbital dynamics. 
The upcoming missions, MESSENGER (Solomon et al. 2001) and BepiCo- 
lombo, with onboard instrumentations capable of measuring the rotational 
parameters of Mercury (see e.g. Milani et al. 2001), stimulate the objective to 
achieve a very accurate theory of Mercury's rotation. One aim of these missions 
is to determine the existence or not of a liquid core by measuring accurately 
both the rotation and the gravity field of Mercury. The current method for 
obtaining constraints on the state and structure of Mercury's core is based 
on the assumptions introduced by Peale in 1972 and revisited in 2002 (Peale 
et al. 2002). By measuring the amplitude of the 88-day libration in longitude, 
called (p hereafter, the obliquity f] (a constant derived from the Cassini state 
assumption for Mercury's rotation and playing the role of a mean obliquity), 
and the degree-two gravitational coefficients (20 and C22), it is possible to 
determine the radius of the expected fluid core. 

Let A, B, and C be the principal moments of inertia of Mercury, C m the 
mantle's moment of inertia, M the mass of Mercury, and R its equatorial 
radius. The information on the core size is inferred from the following relation 
(Peale et al. 2002) : 

\B-A)\ MR 2 ) \ C ) C ~ K } 

where a ratio of C m /C equal to 1 should indicate a core firmly coupled to the 
mantle and most likely solid. In the other hand, if the interior of Mercury is 
partly fluid, the ratio is smaller than 1. 

The objective of the rotation parameter measurements consists in restraining 
each factor of the expression (1). The first factor is inversely proportional to 
the (p rotation angle (Peale et al. 2002). The second factor in (1) is equal to 
4C22 while the third factor is related to fj by the following relation (Wu et al. 
1995) : 

MR 2 _ fi sin (I + fj) 

C n sin 7/ [(1 + cos 7/) G201C22 - (cost]) G210C20] 

Rambaux) . 
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where n is the mean orbital motion of Mercury and /i a constant of uniform 
precession of the orbit pole. G 2 oi an d G 2 io are the eccentricity functions of 
Kaula (Kaula, 1966). G201 = | e — re P resen t the first two terms in a series 
expansion where e is the eccentricity while G 2 io = (1 —e 2 )~i. The formula (2) 
expresses the Cassini state of Mercury (Colombo 1966; Peale 1969) where / is 
the inclination of the orbit pole of Mercury relative to the Laplace plane pole 
around which the orbit pole precesses at an uniform rate. Let us recall that in 
a Cassini state, the orbital and rotational parameters are indeed matched in 
such a way that the spin pole, the orbit pole, and the Laplacian pole remain 
coplanar while the spin and orbital poles on average precess at the same rate. 

Table 1 is here 

The rotational motion of Mercury is mainly characterized by a 3:2 spin-orbit 
resonance and a spin-orbit synchronism where axes of precession cones, spin 
and orbital, on average precess at the same rate. In the 2-body problem (Sun 
and Mercury), by using the 3:2 spin-orbit resonance and a few corollary ap- 
proximations in the classic Euler equations of the 3-D rotational motion of 
Mercury disturbed by the solar torque, one may get two characteristic fre- 
quencies, namely the one in longitude $, which is solution of the third Eule- 
rian equation and the other in latitude of the body which is solution of the 
two first Eulerian equations. A Hamiltonian approach of the rotational mo- 
tion of Mercury has been expanded by D'Hoedt & Lemaitre (2004). The two 
resulting frequencies $ and usually called proper frequencies, calculated by 
this analytical approach and those obtained with our SONYR model are in a 
good agreement (see Table 1). By the present paper, we have also improved 
the internal accuracy of our values (see Table 1). 

Table 2 is here 

Let us call 77 the instantaneous obliquity of Mercury, defined as the angle bet- 
ween Mercury's Oz spin axis and the normal to its orbital plane; let < 77 > 
be its mean value. Until the Radar observations recently obtained by Margot 
et al. (2006) Qthe value of Mercury's mean obliquity was not truly constrai- 
ned by measurements; many values were available in literature from zero to 7 
arcminutes (amin hereafter). On the one hand, from the formula (2) computed 
with realistic Mercury's parameters, fj can not be equal to zero. On the other 
hand, 7 amin given by Wu et al. (1995) seemed too high for Mercury. In our 
previous paper (Rambaux & Bois 2004), using (2) we have estimated a range 
of fj from 1.33 to 2.65 amin (for / 6 [5, 10] deg, e G [0.11,0.24], with [i = 2n/U 
and n = 2ir/P\ where values of the II orbital precession period and the P\ 
orbital period are given in Table 2). However by using our SONYR model 
including the full spin-orbit motion of Mercury, for the plausible values of 

1 published at the very time when we are finishing the revised paper. 
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parameters we had chosen, the mean value of r] is close to 1.66 amin. In the 
present paper, by improving our method, we give a more relevant value for 
< 77 > that proves to be in a very good agreement with the observations. 
We obtain the dynamical behavior of the hermean obliquity by the way of the 
following relation : 

cos 77 = cos i cos 9 + sin i sin 9 cos(f2 — ifj) (3) 

where the spin-orbit variables are computed with SONYR. % and Q are respec- 
tively the inclination and the longitude of the ascending node of Mercury's 
orbital plane relative to the ecliptic plane. The angles ip and 9 are the preces- 
sion and nutation angles of the Eulerian sequence 3-1-3, while the third angle 
ip is called the proper rotation. This angular sequence is used to describe the 
evolution of the body-fixed axis system Oxyz (centered on the Mercury's cen- 
ter of mass) relative to a local reference system, dynamically non-rotating but 
locally transported with the translational motion of the body. Let OX{Y\Z\ be 
the local reference system for the Mercury's rotation (see complete definitions 
of reference systems in Section 3). 

In our previous work as well as in the present one, using our SONYR model, 
we have numerically integrated the complete equations of the iV-body problem 
(8 planets and the Sun) including the full spin-orbit motion of Mercury. More 
precisely, in our previous paper (Rambaux & Bois 2004), we presented (i) an 
accurate theory of Mercury's spin-orbit motion, (ii) the proper frequencies 
related to the spin-orbit resonance, (iii) an analysis of main librations, and 
(iv) a first attempt for characterizing the dynamical behavior of Mercury's 
obliquity. The present paper is wholly devoted to stress and highlight the 
fourth point. 

For lack of current accurate observations or ephemerides of Mercury's rotation, 
and therefore for lack of valid initial conditions for a numerical integration, we 
have built an original method for finding the libration center of the spin-orbit 
system and as a consequence for avoiding arbitrary amplitudes in librations of 
the spin-orbit motion as well as in Mercury's obliquity. Initial conditions of the 
rotational and orbital modes have indeed to be consistent between themselves 
with respect to a same spin-orbit reference system. As a consequence, using 
the SONYR model, the method requires in a first step to build mean initial 
conditions verifying geometrical conditions of a Cassini state for Mercury, 
that is to say a spin-orbit equilibrium state. In a second step, these mean 
initial conditions are fitted in order to place the spin-orbit system at its center 
of libration. The method has been carried out in two cases: (1) the spin- 
orbit motion of Mercury in the 2-body problem case (Sun-Mercury) where an 
uniform precession of the Keplerian orbital plane is kinematically added at a 
fixed inclination, (2) the spin-orbit motion of Mercury in the iV-body problem 
(Sun and planets). With the spin-orbit system located at its center of libration 
in both cases above, the reference system of the first case is put at the average 
of perturbations of the second case. As a consequence, we make not only in 
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evidence the strict dynamical behavior of Mercury's obliquity in both cases 
but may compare the respective oscillations. 

Our theoretical results are in a good agreement with recent radar observa- 
tions. However, it is not excluded that it should be possible to push farther the 
convergence process by drawing nearer still better to the libration center. In 
the present framework, we note that the dynamically driven spin precession, 
which occurs when the planetary interactions are included, is geometrically 
more complex than the purely kinematic case. In particular, we highlight that 
in the full spin-orbit motion of Mercury, the asynchronous rotation of the 
planet is really combined with a synchronism in precession variables as well 
as a second one between secular changes of the orbit inclination and the spin 
axis orientation. As a consequence and supposing that it could exist an irre- 
ducible minimum of the oscillation amplitude, we show how the quasi-periodic 
oscillations found in Mercury's obliquity could be geometrically understood. 
We discuss briefly the question to know whether the minimal amplitude has 
to theoretically be of zero value or not. In the meantime of such a forthcoming 
work, we describe Mercury's rotation and the resulting equilibrium state that 
proceeds from a 2-body Cassini state. 



2 The Spin-Orbit AT-bodY Relativistic model 



For obtaining the real motion of Mercury, we have used our model of solar 
System integration including the coupled spin-orbit motion of the Moon. This 
model (built by Bois, Journet, & Vokrouhlicky and called BJV), expanded 
in a relativistic framework, had been previously built in accordance with the 
requirements of the Lunar Laser Ranging observational accuracy (Bois 2000; 
Bois & Vokrouhlicky 1995). We extended the BJV model by generalizing the 
Moon's spin-orbit coupling to the other terrestrial planets (Mercury, Venus, 
the Earth, and Mars). The model is at present called SONYR (acronym of 
Spin-Orbit iV-BodY Relativistic model). As a consequence, the SONYR model 
gives an accurate simultaneous integration of the spin-orbit motion of Mer- 
cury. The approach of the BJV and SONYR models consists in integrating the 
A^-body problem (including translational and rotational motions) on the basis 
of the gravitation description given by the general relativity. The equations 
have been developed in the DSX formalism presented in a series of papers by 
Damour, Soffel, & Xu (1991, 1992, 1993). For purposes of celestial mechanics, 
to our knowledge, it is the most suitable formulation of the post-Newtonian 
(PN) theory of motion for a system of N arbitrary extended, weakly self- 
gravitating, rotating and deformable bodies in mutual interactions. The DSX 
formalism, derived from the first PN approximation level, includes then the 
translational motions of the bodies as well as their rotational motions with 
respect to locally transported frames with the bodies. To each body is then 
associated a local reference system, centered on its center of mass and dynam- 
ically non-rotating. These local reference systems are freely falling in the local 
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gravity field (the Sun) and are referred to a global reference system O'X'Y'Z' 
for the complete iV-body problem. 

A description of the SONYR model is given in our previous paper (Rambaux 
& Bois 2004). The principle of the model allows the analysis of Mercury's spin- 
orbit motion and the identification of internal mechanisms such as the direct 
and indirect planetary interactions. In addition, it gives also the interrelations 
between the parameters involved in Mercury's dynamical figure and the planet 
librations. 

3 Reference systems and their relations 

O'X'Y'Z' is centered on the Solar System barycenter, the O'X'Y' plane is 
parallel to the J2000 ecliptic plane. The X'-axis is oriented to the equinox 
J2000. Let OXY be the mean orbital plane of Mercury and OZ the mean 
orbital pole. Let us call V the instantaneous orbital plane of Mercury. The 
rotational motion of Mercury is evaluated from a body-fixed coordinate axis 
system Oxyz centered on Mercury's center of mass relative to a local dy- 
namically non-rotating reference system, OX\Y\Z\. The Ox,Oy,Oz axes are 
defined by the three principal axes of inertia a, b, and c of the planet according 
to the best fitting ellipsoid such as a > b > c. The corresponding principal 
moments of inertia are then ordered as follows : A < B < C. We used the 
Euler angles if), 9, ip related to the 3-1-3 angular sequence (3 represents a ro- 
tation around a Z-axis and 1, around an X-axis) to describe the evolution of 
the body-fixed axes Oxyz with respect to axes of the local reference system 
OX\Y\Z\. These angles are defined as follows : ip is the precession angle of the 
Oz polar axis around the OZ\ reference axis, 9 is the nutation angle represent- 
ing the inclination of Oz with respect to OZ\, and p is the rotation around 
Oz and conventionally understood as the rotation of the largest rotational 
energy. This angle is generally called the proper rotation. The Oz principal 
axis around which is applied the proper rotation is called the axis of figure 
and defines the North pole of the rotation (Bois 1995). The proper rotation <p 
is around the smallest axis of inertia; it is called hereafter the spin angle and 
Oz the spin axis. 

Figure 1 is here 

In SONYR, building up the reference frames from the reference systems is re- 
lated to the choice of the reference frames transported by the initial conditions 
(IC hereafter). This is of practical use when we have valid IC (for instance 
coming from ephemerides or derived from observations). In such conditions, 
the local reference frame for the rotation of a planet may be the one of the 
ephemeris. It is not actually the case for Mercury. True IC of Mercury's ro- 
tation are unknown (with respect to any reference system). Moreover, ad hoc 
attempts for determining instantaneous IC of the rotation produce a shift with 
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those of the orbit and therefore a shift of the libration center of the spin-orbit 
resonance. Inconsistency between rotation's and orbit's IC drives necessarily 
to arbitrary amplitudes and phases in the spin-orbit librations. As a conse- 
quence, we have built (in Section 4) mean initial conditions (MIC hereafter) 
coming from geometrical conditions of the mean spin-orbit motion of Mercury. 
The latter is related to (i) the 3:2 spin-orbit resonance in cyclic variables (if 
and A such as < dtp/dt > = (3n/2) + < du/dt >), (ii) the synchronism in 
precession variables (ip and Q), and (iii) a second synchronism in nutation 
and inclination variables (9 and i) (see Rambaux & Bois 2004). In a second 
step of the method (Section 5), the MIC are accurately corrected by taking 
into account a fit between the central axis of the precession cone described by 
the spin pole and the theoretical axis of reference. For these final MIC, the 
spin-orbit system is at its center of libration (r hereafter) without arbitrary 
amplitudes. Moreover, we consider two cases : 

(1) The spin-orbit motion of Mercury in the 2-body problem framework (Sun- 
Mercury) where a uniform precession, which we denote by II, of the Keplerian 
orbital plane is kinematically added at a fixed inclination Ik with respect to 
OX'Y' . Let us call S2K this construction bearing the geometrical Cassini state. 
Let us call Vk the associated orbital plane (a Keplerian plane in this case), 
and Ck the fixed plane of reference OX'Y' (Ck plays the role of a Laplace 
plane for Vk) ( see Fig- !)• 

(2) The spin-orbit motion of Mercury in the iV-body problem framework (Sun 
and planets). Let us call S n this case. Let us call V n the associated mean orbital 
plane of Mercury. The instantaneous orbital plane V has an instantaneous 
inclination % relative to OX'Y' while the inclination of V n is < i >. Let us 
include the S2K construction into S n by putting Vk on V n . As a consequence, 
Ik = < i >• Over short periods, V n remains equivalent to Vk- The associated 
Laplace plane C n is an instantaneous plane which follows the oscillations of V 
so that i remains quasi- const ant relative to C n . Ck coincides with the average 
of the oscillations of C n , i.e. < C n > = Ck (see Fig. 1). 

Our construction of reference systems verifies in the end the three following 
relations : 

V n = V K <t>=I K <C n >=C K (4) 

which permit to accurately evaluate the obliquity's oscillations (i) not only 
with respect to the center of libration of the full spin-orbit motion of Mer- 
cury, but (ii) with respect to the included Cassini state. In other words, our 
construction is such as the mean Cassini state of the full spin-orbit motion 
of Mercury, on a given interval of time, coincides with an exact Cassini state 
derived from the 2-body problem. 

Let us add that local dynamically non-rotating frames show in general relativ- 
ity a slow (de Sitter) rotation with respect to the kinematically non-rotating 
frames. As a consequence, the OX{Y\Z\ local reference frame for Mercury's 
rotation is affected by a slow precession of its axes transported with the trans- 
lational motion of Mercury. In the Earth's case, the de Sitter secular preces- 
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sion of the Earth reference frame is close to 1.9"/cy (see Fukushima 1991; 
Bizouard et al. 1992; Bois & Vokrouhlicky 1995). Consequently, the real rota- 
tion of Mercury has not to be expressed in an inertial system fixed in space, 
but in a local dynamically non-rotating frame fallen down in the gravitational 
field of the Sun. However, in the framework of the present paper where only 
Newtonian-like terms of SONYR are activated, OX\Y\Zi remains parallel to 
OX'Y'Z'. 



4 Building the mean initial conditions 



For lack of accurate observations or ephemerides of Mercury's rotation, and 
therefore for lack of valid IC for a numerical integration of the full spin-orbit 
motion of Mercury, a method is required for: (i) fitting the different parameters 
or constants involved in the spin-orbit of Mercury so that the SONYR model 
be self-consistent; (ii) removing arbitrary amplitudes in librations of Mercury's 
spin-orbit motion and therefore in its obliquity. In this section, we present the 
construction of MIC related to the spin-orbit coupling of Mercury. 

(1) We start from the geometry of the spin-orbit coupling problem defined in 
Goldreich & Peale (1966): the spin axis of Mercury is considered normal to 
its orbital plane and the orbit is assumed to be fixed and unvarying. At this 
level, for expressing these assumptions, as well as that the long axis is pointed 
toward the Sun at a perihelion position to, analytical relations giving the 
corresponding set of initial conditions, namely Co, can be written as follows : 



tp = u + f M 



(C 






where ujq is the mean argument of periastron, and M the mean anomaly with 
Q = loq = M — at t — to- For such a Mercury's spin-orbit coupling based 
on a 2-body problem (with the asphericity factor a = \J3(B — A)/C = 0.0187 
and the eccentricity e = 0.206), we have plotted a Poincare section in Rambaux 
& Bois (2004, Fig. 2), making in evidence the quasi-periodic librations around 
the 3:2 rotation state. 

(2) In the S2K framework, the condition of the precession synchronism drives 
to the geometrical property < dip/dt >= 2tt/H where II is the constant of 
uniform precession. The principle of our method requires to obtain the value of 
II from the S n case such as 2tt/U = < dQ/dt >. In the same way, the value of 
I = Ik is obtained from the S n case such as Ik = < i >, and as a consequence 
the value of / is taken equal to < i > since the S2K case. The mean initial 
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conditions Cik for the Sik case are then written as follows 



ip = tt 



¥o = u + |M 



(C2K) 




2tt 

TT 





with f2 , to , and M equal to zero at t 

(3) In the S n framework, the orbit precession of the rotating body is dy- 
namically generated by the planetary interactions but its rate is nonuniform. 
However, 2n/U has been rightly defined as < dVt/dt >. Moreover, the 3:2 spin- 
orbit resonance in the S n case involves the condition dcp/dt = (2n/A) + (3n/2) 
where A is obtained such as 2ir/A = < duj/dt >. In addition, as shown further, 
9 and i on average change at the same rate. This expresses a proximity to a 
Cassini state, which involves a third geometrical condition for building MIC, 
namely < dO/dt > = < di/dt >= 2n/Q. The mean initial conditions C n for 
the S„ case are then written as follows : 



(Cn) 



V> = ^0 



6c 




2tx 




< % > 

2tx 



(p =UJ + 



1M 



d(p 



2ix 3 
= — — I — n 
dt J Q A 2 



with fl , lo , and M equal to zero at t = to. 

The dynamical behaviors of 9 and i have the same slopes (see Fig. 2, top 
panel). As a consequence, secular parts of 9 and i decrease at the same rate 
equal to 6. Although it is a question of a A^-body problem, the spin keeps its 
proximity to a Cassini state during the slow change in i. But with this signature 
of the proximity to the Cassini state made in evidence in the S n framework, 
< d9/dt > = < di/dt > may be understood as a second synchronism between 
the respective inclinations of the orbit and the spin axis relative to a same 
reference system (OX'Y' Z'). As a matter of fact we find that £ = 9 — i librates 
and rightly according to the \l/ proper frequency in latitude (see Fig. 2, bottom 
panel) . 



Figure 2 is here 



5 Fitting the mean initial conditions 

Let us call hereafter SONYR(5*j, Cj) the integration of SONYR in the mode 
Si with the initial conditions CV Let us note that under the assumptions of 
the present paper, SONYR(S'j, Ci) is a conservative dynamical system. Let S 
be the Oz spin pole of Mercury. S describes a large cone with oscillations 
of polhodie; the projection of the cone in the OX'Y' reference plane gives a 
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large quasi-circle of center A (with oscillations). Due to slightly heterogeneous 
initial spin-orbit positions, A does not exactly coincide with the projection of 
OZ' . By setting the latter in A, the system is placed at its V libration center. 
This can be achieved by slightly moving OX'Y'Z' or in SONYR by fitting a 
little bit the MIC. But because of large periods of precession, building and 
fitting the large circles is a very heavy computational process. Nevertheless, 
the reference frame that satisfies the whole system in V is the one where all 
librations (related to small and large cones of precession) are at their respective 
right centers of libration. As a consequence, we have determined equations 
expressing the departure from V within smaller circles of libration. Let us note 
Sx" and Sy the coordinates of S projected onto a reference plane OX"Y" 
rotating with the uniform motions IT and 0; they are written as follows : 

S x »= sin -(£)*] sin [0 -(§)*] 
Sy „ = _ CO s[^-(2E)*]sm[0-(f)t] 

where 2n/0 = in the S 2 k case. For each mode of SONYR (S 2 k and S n ), 
S describes a small precession cone. Figure 3 gives the motion of S in the 
mode SONYR(S , 2x, C 2 k). Let us call 0( the central axis of the small cone of 
precession whose coordinates are (Ci,Ca) i n OX"Y" . 



Figure 3 is here 



The fitting equations are then written as follows : 

Sx - = ° (6) 
Sy — C2 

After substitutions in orbital elements alone, the system (6) becomes : 



Ik + & Ik = arcsin (y Ci + Cl) ^ 
Q + 5Q = arctan (— C1/C2) 

where one gets SIk and SQ as functions of Ik, fio ; Ci> and ( 2 - In the mode 
SONYR(S , 2_ft-, C 2 k), we obtain then the two following fitted initial conditions : 

ipo = tt + 5tt 

e = i K + 51 K 

C 2 k becomes C' 2K and one gets on the same graphic frame as previously a 
little circle (like a point at this scale, see Fig. 4). The little circle expresses 
a considerable decreasing of the degree of opening of the cone by the above 
process. 



Figure 4 is here 
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The same principle is applied to the S n case with an analogous result. Taking 
into account in C n (becoming C' n ) the two fitted initial conditions ipo and 8q 
coming from (8) according to the same process, Figure 5 shows the significant 
decreasing of the degree of opening of the precession cone. 

Figure 5 is here 

Table 3 is here 



6 The resulting obliquities 

In the S 2 k mode, computations of r] give resulting obliquities before and after 
fitting the MIC at the libration center T, as shown in Figure 6a. Figure 6b 
shows a zoom of the obliquity purged of arbitrary amplitudes after the fit- 
ting process. However, rj appears not constant but quasi- const ant. In the S 2 k 
mode indeed, i.e. without planetary interactions, due to the gravitational so- 
lar torque acting on the dynamical figure of Mercury with three unequal axes 
of inertia, the equations of the rotational motion solely form a non-integrable 
system whose regular solutions are quasi-periodic librations. As a consequence, 
the rj obliquity that depends on 9 and ip (see formula (3) while i = Ik and 
Q = Q + (27r/II)t where Ik, an d II are constant values) follows in part 
this behavior. Nevertheless, with such faint amplitudes (peak-to-peak 0.4 arc- 
seconds, 'as' hereafter), the obliquity may be reported quasi- const ant in the 
S 2 k case, which is usually characteristic of a Cassini state derived from the 
2-body problem (see Fig. 6b). After removing the amplitude of the high fre- 
quency fluctuations (a modulation of a beat frequency with a period of 295 
years whose peak-to-peak amplitude varies between a maximum near 0.28 as 
and a minimum near 0.14 as), the peak-to-peak amplitude of the resulting 
long period modulation of 1066.4 years is only 0.12 as. 

Figure 6 is here 

In the same way, in the S n mode, computations of rj give resulting obliquities, 
before and after the fitting process, as shown in Figure 7a. Figure 7b shows 
a zoom of the obliquity after the fitting process and then purged of arbitrary 
amplitudes. In both modes, S 2 k and S n , the period of the main oscillation is 
namely the second proper frequency of the resonance. A frequency analysis 
gives 1066.81 years in the S n case and 1066.44 years in the S 2 k one, let be a 
difference lesser than 0.4 years. In the S n case, after removing the amplitude 
of the high frequency fluctuations (about 0.6 as), the peak-to-peak amplitude 
of the resulting long period modulation of 1066.8 years is 3.4 as. 

Figure 7 is here 
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Both sets of results, from S 2 k and S n cases are globally consistent with those 
of Peale (2006) and implicitly with those of Yseboodt & Margot (2006). We 
note that the dynamically driven spin precession, which occurs when the plan- 
etary interactions are included, is geometrically more complex than the purely 
kinematic case. We find that the 3:2 spin-orbit resonance is not only combined 
with a first synchronism where the spin and orbit poles on average precess at 
the same rate but also with a second one where the rotational nutation desi- 
gnated by 9 and the i orbit inclination on average decrease at the same rate 
(see Fig. 2 where is shown the libration of £ = 9 — i within the \1/ period). 
These two synchronisms reflect the proximity to the Cassini state and their 
signatures express its conserving in the TV-body problem case. Besides, the 
conceptual question to know whether the ultimate minimal amplitude could 
be of zero value or not is discussed further. However, regardless of that, inside 
the present numerical process, the amplitudes in the S n case are related to 
a spin-orbit coupling. As shown below with the equation (9), it is about an 
indirect effect due to the solar torque, bearing the planetary interactions, and 
acting on the planet's rotation. 

Reference frames of the S 2 r an d S n cases being linked by the conditions (4) 
(see Fig. 1) and values of the free constants required for SONYR(S , 2 x, C 2 k) 
being determined from SONYR(5 n , C n ), SONYR(SW, C 2 k) has been built 
at the heart of SONYR(S' n , C n ) in such a way the quasi-periodic solution 
of SONYR(S' n , C n ) surrounds the vicinity of the quasi-periodic solution of 
SONYR(S < 2_ft-, C<ik)- Arbitrary amplitudes are not only avoided in the oscilla- 
tions of rj, but a direct comparison of 7] coming from S n with r] coming from S 2 k 
is possible, which makes in evidence the effect of the planetary interactions. 
In the S n mode indeed, the gravitational solar torque, bearing the planetary 
interactions and acting on the dynamical figure of Mercury, induces more size- 
able librations related to the disturbed Euler angles as well as the variables 
i and Q. At this phase of analysis, the obliquity can be no longer said quasi- 
constant. The oscillation amplitude of period \1/ is indeed equal to 4 as (total 
peak-to-peak amplitude), let be one order of magnitude larger than in the 
S2K case (or 28 times larger when comparing the peak-to-peak amplitudes 
of the strict \P modulation without high frequency fluctuations) (see Fig. 6b 
and Fig. 7b). The comparison of the S 2 k and S n modes, respectively derived 
from the 2-body and iV-body problems, both computed and plotted at their 
respective libration centers T, shows how Mercury's obliquity and its dynami- 
cal behavior are modified when the planetary interactions are fully introduced 
(see Fig. 8). The mean obliquity is larger as well as the oscillation amplitude. 
< 77 > = 1.98 amin in the S n case while < 77 > = 1.80 amin in the S 2 k case. 

This spin-orbit mechanism is rightly an indirect effect due to the solar torque, 
bearing the orbital planetary interactions and acting on Mercury's rotation, 
as shown by the classic general equation of the rotational motion of a solid 
celestial body (Bois 1995) : 
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<9[(i>] /3Gm\ /a\ 3 

-^p + a; x [(/M = (-) u x (/) tt (9) 

where oj is the instantaneous rotational vector of Mercury and (J) its tensor 
of inertia. The right hand side of the equation represents the second degree 
torque exerted by a point mass m (the Sun) of position ru with respect to 
the body (Mercury) when it is reduced to three oblateness coefficients, a is 
the mean distance between the two bodies while r is the instantaneous one. 
ru is solution of the iV-body problem in SONYR(S' n ). Let us note that the 
direct planetary torques are about two orders of magnitude smaller than the 
indirect effects acting by the solar torque. 



Figure 8 is here 



Due to the planetary interactions, Mercury's orbit is no longer Keplerian and 
rate of precession of Mercury's orbital plane is no longer uniform. As a conse- 
quence, the obliquity's dynamical behavior is no longer quasi-constant. More 
generally, it is not possible for the obliquities to be constant since the orbit 
precession occurs at non-uniform rates (Bills 2005). In the case of Mercury, 
according to Bills & Comstock (2005), if the orbit pole precession rate were 
uniform, dissipation would have driven Mercury into a Cassini state, in which 
the spin pole and orbit pole remain coplanar with the Laplace pole, as the 
spin pole precesses about the moving orbit pole. However, according to these 
authors and in accordance with our results, due to the nonuniform orbit pre- 
cession rate with several distinct modes of oscillation, this simple coplanar 
configuration is no longer globally maintained. 



7 Discussion 



We have significantly reduced amplitudes of the oscillations in Mercury's obli- 
quity, which were exposed to discussion since a previous paper (Rambaux & 
Bois 2004). For doing that, we have built an original method for finding the 
libration center of the spin-orbit system and as a consequence for avoiding ar- 
bitrary amplitudes in librations of the spin-orbit motion as well as in Mercury's 
obliquity. We have converged and the amplitudes are in very good agreement 
with most recent observations (Radar observations by Margot et al. 2006). 
This validates not only our theoretical SONYR model applied to Mercury but 
also our method of convergence. 

Our centering method has converged but it is not excluded that it may re- 
main a residual amplitude of arbitrary oscillations. It ought to be possible to 
push farther the convergence process by drawing nearer yet more precisely 
to the libration center (in the framework of the paper, i.e. in a conservative 
dynamical system). Since our amplitudes are compatible with current observa- 
tions, we have concluded that this question could wait for MESSENGER and 
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BepiColombo data in order to directly fit the SONYR model to observations. 
Nevertheless, apart from values of residual amplitudes, it remains a conceptual 
question. Does there exist a theoretical possibility to reach a pure point of 
equilibrium where the oscillation amplitude would be absolutely of zero value 
or, on the contrary, does there exist an irreducible minimum, even very small 
but of non-zero value? In the S n case, existence of distant resonances of high 
orders that would excite the 3:2 resonance is not excluded a priori. None 
mathematical theorem does deny it but the demonstration of this would take 
us along to a purely mathematical ground. Consequently, we will say that the 
question remains open (both in conservative and dissipative cases). However, 
we plane some further works in this connection. The occurrence of a perfect or 
rigid equilibrium without oscillations seems to us unlikely for natural systems. 

8 Conclusion 

Using our SONYR model of the solar System including the spin-orbit cou- 
plings of Earth-like planets, and setting the system at its libration center in 
order to eliminate arbitrary amplitudes (at the internal accuracy level of the 
model), we have obtained the complete obliquity for the full spin-orbit mo- 
tion of Mercury. Our construction of reference systems permits to accurately 
evaluate the obliquity's oscillations (i) not only with respect to the libration 
center of the full spin-orbit motion of Mercury (S n ), (ii) but also with respect 
to a Cassini state derived from the 2-body problem (SW)- The construction 
is such as the exact Cassini state derived from S2K coincides with the mean 
Cassini state of S n on a given interval of time. We have obtained (1) a mean 
value of Mercury's obliquity that is in a very good agreement with very recent 
observations (Radar observations by Margot et al. 2006), namely 1.98 amin 
with SONYR for 2.1 ±0.1 amin by observations, (2) amplitudes of oscillations 
in Mercury's obliquity of the order of 2 as, which is compatible with the results 
of these observations where the deviation with respect to the 2-body Cassini 
state is constrained lesser than 6 as. These observational data validate not 
only our theoretical prediction with the SONYR model applied to Mercury 
but also our method of convergence for finding a libration center in a complex 
spin-orbit system. 

Being dependent on the question of existence of an irreducible minimum in the 
oscillation amplitude, we find that the equilibrium state of Mercury computed 
in the framework of a conservative iV-body system squares to a quasi-periodic 
oscillation without secular drifts. In other words, in spite of the planetary 
interactions indirectly acting by the solar torque on Mercury's rotation, Mer- 
cury remains in a stable equilibrium state that proceeds from a 2-body Cassini 
state. 

Besides, by the way of our method, we have made in evidence the richness of a 
dynamical precession related to planetary interactions (iV-body problem) rel- 
ative to a purely kinematical orbit precession added to a 2-body problem (Sun 
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and Mercury). Whatever the open question on the minimal amplitude in the 
obliquity's oscillations, the approach of the 2-body Cassini state is necessarily 
limited to reach a constant or quasi- const ant behavior of Mercury's obliquity 
(where moreover < 77 > = 1.80 amin versus 1.98 amin, for a difference of 
10.8 as). The planetary interactions generate a forced dynamical precession 
of Mercury's orbital plane with a nonuniform rate. This precession is kine- 
matically taken into account in the conventional construction of a Mercury's 
Cassini state as well as in our S2K case. But in the S n case, due to the first 
synchronism where the spin and orbit poles on average precess at the same 
rate (as well as the second one where the spin axis orientation, namely the 
9 nutation, and the orbit inclination % on average change at the same rate), 
the difference angle between these two poles, i.e. rj, (and 9 — i as well) li- 
brates necessarily. As a consequence and whether it would turn out that there 
exists an irreducible minimum of the oscillation amplitude, so these resulting 
quasi-periodic oscillations in the obliquity ought to be understood as librations 
related to these two synchronisms where the (ip — Q) and (9 — 1) respective 
differences librate (according to the \l/ proper frequency in latitude). In the 
same way, in multi-planet systems, when the longitudes of periapse on average 
precess at the same rate, the difference angle between the two apsidal lines 
librates (see e.g. Lee & Peale 2003 or Bois et al. 2003). Under this assumption 
and according to the obtained dynamical behavior, the obliquity of Mercury 
should be characterized by a quasi-periodic oscillation according to the period 
with an amplitude A v and a mean obliquity < 77 >, in such a way that 
Mercury's obliquity 77 could be written as follows (in first approximation) : 



where a is a phase factor. 

In the end, for lack of current and accurate observations or ephemerides of 
Mercury's rotation, but with SONYR(5' ri , C' n ) set up at the libration center 
of Mercury's spin-orbit motion, we may conclude that our method is relevant 
to link physical causes and their respective effects. As a consequence, with 
accurate observations of Mercury's rotational motion, that should be being 
achieved by MESSENGER and BepiColombo, SONYR^, C' n ) , and not a 2- 
body Cassini state model as SONYR(5 , 2^, C' 2K ), ought to be fitted to these 
observations. 
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Proper Frequencies 



$ (years) 



^ (years) 



Rambaux k Bois 2004 15.847 1066 

D'Hoedt & Lemaitre 2004 15.857 1066.68 

This paper 15.824 ± 0.024 1066.8 ± 0.4 

Table 1 

Proper frequencies related to the spin-orbit motion of Mercury (in longitude: <3?, in 
latitude: ty). Comparaison of results. 
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Basic periods 


P v (rotation period) 


= 58.646 days 


Pa (orbital period) 


= 87.969 days 


II (orbital precession) 


= 279 000 years 



Table 2 

Basic periods of the spin-orbit motion of Mercury (Rambaux & Bois 2004) 
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S2K 


S n 


5Qq (amin) 


0.0 


-6.564 


51k (amin) 


-1.806 


-1.823 



Table 3 

Numerical values of 5£Iq and 51k in the S2K and S n cases. 
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Fig. 2. Second synchronism of the spin-orbit motion of Mercury. 9 and i are plotted 
in the top panel over 1600 years. The difference £ = 9 — i plotted on the bottom 
panel librates within the proper frequency in latitude, namely ^ = 1066.8 years. 
Degrees (top panel) and arcminutes (bottom panel) are on the respective vertical 
axes. Years are on both horizontal axes. 
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Fig. 4. Projection of Mercury's spin pole S on the OX"Y v plane. The large 
circle (1) obtained in mode SONYR^S^, Cik) is the one of Fig. 3. In mode 
SONYR(S , 2_ftr, C' 2K ), it is a broadly reduced circle (2). 
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Fig. 5. Projection of Mercury's spin pole S on the OX"Y" plane. The large circle 
(1) is obtained in mode SONYR(S' n , C n ) whereas the broadly reduced circle (2) in 
mode SONYR(S n ,C;). 
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Fig. 6. (a) Resulting obliquities by SONYR(5 2 ^, C 2K ) and SONYR(S 2 k, C' 2K ). 
(b) Zoom of rj in the SONYR(S , 2x, C 2K ) case, i.e. without arbitrary amplitudes. r\ 
is characterized by a main oscillation according to the second proper frequency \& 
within a peak-to-peak amplitude lesser than 0.4 as. In both panels, arcminutes are 
on the vertical axes and years on the horizontal ones. 



28 



2.01 



200 400 600 800 1000 1200 1400 1600 
T (years) 



(b) 














































































1 




i 





1.99 
1.98 
1.97 
1.96 
1.95 
1.94 
1.93 

200 400 600 800 1000 1200 1400 1600 
T (years) 

Fig. 7. (a) Resulting obliquities by SONYR(5 n , C n ) and SONYR(,S n , C' n ). (b) Zoom 
of rj in the SONYR(S' n , C' n ) case, i.e. without arbitrary amplitudes, r/ is charac- 
terized by a clear oscillation according to the second proper frequency while the 
peak-to-peak amplitude is near 4 as. Arcminutes are on the vertical axes and years 
on the horizontal ones. 
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Fig. 8. Comparison of the obliquity's dynamical behavior in SONYR(S' n , C' n ) and 
SONYR(52_ftr, C 2K ) modes, i.e. with and without planetary interactions, both com- 
puted and plotted at their respective libration center V. < r\ > = 1.98 amin in the 
S n case while < rj > = 1.80 amin in the S2K case. 
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